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Abstract 

We study a porous medium equation with fractional potential pressure: 

d t u = V-(u m ~ 1 Vp), p = (-A)~ s u, 

for m > 1, 0 < s < 1 and u(x,t) > 0. The problem is posed for x £ R w , N > 1, and t > 0. 
The initial data u(x, 0) is assumed to be a bounded function with compact support or fast decay 
at infinity. We establish existence of a class of weak solutions for which we determine whether the 
property of compact support is conserved in time depending on the parameter to, starting from the 
result of finite propagation known for m = 2. We find that when to £ [1, 2) the problem has infinite 
speed of propagation, while for m £ [2, 3) it has finite speed of propagation. In other words m = 2 
is critical exponent regarding propagation. The main results have been announced in the note |29| . 
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1 Introduction 


In this paper we study the following nonlocal evolution equation 


( 1 . 1 ) 


u t (x,t) = V • ( u m 1 Vp), p = (—A) s u, for x £ R N , t > 0, 
u( 0, x) = uq(x) for x £ R N , 


for to > 1 and u(x, t ) > 0. The model formally resembles the classical Porous Medium Equation (PME) 
Ut = A u m = V(rmt m_1 Vu) where the pressure p depends linearly on the density function u according to the 
Darcy Law. In this model the pressure p takes into consideration nonlocal effects through the Inverse Fractional 
Laplacian operator K, s = (—A) _s , that is the Riesz potential of order 2s. The problem is posed for x £ R N , 
N > 1 and t > 0. The initial data uq ■ R N —> [0, oo) is bounded with compact support or fast decay at infinity. 

As a motivating precedent, in the work m Caffarelli and Vazquez proposed the following model of porous 
medium equation with nonlocal diffusion effects 


(CV) 


d t u = V • (uVp), p = {— A) s u. 


The study of this model has been performed in a series of papers as follows. In [10], Caffarelli and Vazquez 
developed the theory of existence of bounded weak solutions that propagate with finite speed. In m . the same 
authors proved the asymptotic time behaviour of the solutions. Self-similar non-negative solutions are obtained 
by solving an elliptic obstacle problem with fractional Laplacian for the pair pressure-density, called obstacle 
Barenblatt solutions. Finally, in [Sj, Caffarelli, Soria and Vazquez considered the regularity and the L 1 — L°° 
smoothing effect. The regularity for s = 1/2 has been recently done in !) . The study of fine asymptotic 
behaviour (rates of convergence) for ( |CV| ) has been performed by Carrillo, Huang, Santos and Vazquez )T2] in 
the one dimensional setting. Putting to = 2 in (|1.1[), we recover Problem (|CV|). 


A main question in this kind of nonlocal nonlinear diffusion models is to decide whether compactly supported 
data produce compactly supported solutions, a property known as finite speed of propagation. Surprisingly, the 
answer was proved to be positive for to = 2 in paper m, for m = 1 we get the linear fractional heat equation, 
that is explicitly solvable by convolution with a positive kernel, hence it has infinite speed of propagation. The 
main motivation of this paper is establishing the alternative finite/infinite speed of propagation for the solutions 
of Problem (1.1) depending on the parameter to. In the process we construct a theory of existence of solutions 
and derive the main properties. A modification of the numerical methods developed in mm pointed to us to 
the possibility of having two different propagation properties. 


Other related models. Equation (CV) with s = 1/2 in dimension N = 1 has been proposed by Head [2D] to 
describe the dynamics of dislocation in crystals. The model is written in the integrated form as 

v t + \v x \(-d 2 /d xx ) 1/2 v = 0. 

The dislocation density is u = v x . This model has been recently studied by Biler, Karch and Monneau in [J], 
where they prove that the problem enjoys the properties of uniqueness and comparison of viscosity solutions. 
The relation between u and v is very interesting and will be used by us in the final sections. 

Another possible generalization of the (CV) model is 

d t u = V-(uVp), p=(-A)- s (\u\ m - 2 u), 

that has been investigated by Biler, Imbert and Karch in HE]- They prove the existence of weak solutions and 
they find explicit self-similar solutions with compact support for all to > 1. The finite speed of propagation for 
every weak solution has been done in [122] , 

The second nonlocal version of the classical PME is the model 

m = -(-A) s 'it m ', to'>0, 
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known as the Fractional Porous Medium Equation (FPME). This model has infinite speed of propagation and the 
existence of fundamental solutions of self-similar type or Barenblatt solutions is known for m > ( N — 2 s') + /N. 
We refer to the recent works on cm on Eg. The (FPME) model for nn! = 1, also called linear fractional Heat 
Equation, coincides with model for s = 1 — s', m = 1. 


1.1 Main results 

We first propose a definition of solution and establish the existence and main properties of the solutions. 


Definition 1.1. Let m > 1. We say that u is a weak solution of (1.1) in Qt = R N x (0,T) with nonnegative 
initial data Uq £ L 1 (K Ar ) if (i) u € L 1 (Qt), (ii) V/C s [u] £ L 1 ([0,T ) : L( oc (R n )), (Hi) u m ~ 1 S7IC s [u] £ L 1 (Qx), 
and (iv) 


( 1 . 2 ) 


f j ucftdxdt— f ( u m 1 S7K, s {u)S7(j)dxdt + f uo(x)<j>(x, 0)dx = 0 
J o J i" Jo J r jv Jr n 


holds for every test function </> in Qt such that X/cj) is continuous, (f> has compact support in R N for all t £ (0, T ) 
and vanishes near t = T. 

Before entering the discussion of finite versus infinite propagation, we study the question of existence. We have 
the following result for 1 < m < 2. 


Theorem 1.2. Let m £ (1,2), N > 1. Let uq £ L 1 (R n ) ni°°(] R N ). Then there exists a weak solution u of 
equation (1.1) with initial data uq such that u £ L 1 {Qt) D L°°(Qt) and X7TL s [u\ £ L 2 (Qt)- Moreover, u has the 
following properties: 

1. (Conservation of mass) For all t > 0 we have / u(x,t)dx = / uo{x)dx. 

J r n Jr n 


2. (L°° estimate) For all t > 0 we have ||u(-,t)|| 00 < 

3. (First Energy estimate) For all t > 0, 


cf / \X7Hs[u}\ 2 dxdt+ [ u(tf- m dx< 

Jo Jr n J r w 


ul~ m dx, 


with C =l (2 m) (3 m) > 0. 

4- (Second Energy estimate) For all t > 0, 




\H s [u]\ 2 dx + / [ u m 1 \X7IC s [u]\ 2 dx < ^ f \TL s [u 0 }\ 2 dx. 
Jo JR N ^ Jr n 


The existence for m > 2 is covered in the following result. 

Theorem 1.3. Let m £ [2,3). Let Uq £ L 1 (R N ) fl L°°(R N ) be such that 
(1.3) 0 < uo(x) < Ae~ a \ x \ for some A, a > 0. 

Then there exists a weak solution u of equation with initial data uo such that u £ L 1 {Qt ) D L°°(Qt), 
Vlf s [rt] £ L 2 (Qt ) and u satisfies the properties 1,2,4 of Theorem 1.2. Moreover, the solution decays exponen¬ 
tially in |a:| and the first energy estimate holds in the form 


\C\ [ [ \XUs[u]\ 2 dxdt 

Jo Js. N 

where C = C{m) = (2 — to)(3 — to). 


+ 


,3 —m 


dx < 


/R» 


i(tf~ m dx 


4 





We should have covered existence in the whole range m > 2 where we want to prove finite speed of propagation 
for the constructed weak solutions, see Theorem |1.4| But the existence theory used in the previous theorem 
breaks down because of the negative exponents 3 — m that would appear in the first energy estimate for m > 3 
(a logarithm would appear for m = 3). A new existence approach avoiding such estimate is needed, and this 
can be done but is not immediate. We have refrained from presenting such a study here because it would divert 
us too much from the main interest. 

The following is our most important contribution, which deals with the property of finite propagation of the 
solutions depending on the value of m. 


Theorem 1.4. a) Let N > 1, m £ [2, 3), s £ (0,1) and let u be a constructed weak solution to problem (1.1) as 
in Theorem 1.3 with compactly supported initial data uq £ L 1 (R N ) nL°°(R N ). Then, u(-,t) is also compactly 
supported for any t > 0, i. e. the solution has finite speed of propagation. 

b) Let N = 1, m £ (1,2), s £ (0,1) and let u be a constructed solution as in Theorem 1.2 Then for any t > 0 
and any R > 0, the set M. R t = {x : |x| > R , u(x,t) > 0} has positive measure even ifuo is compactly supported. 
This is a weak form of infinite speed of propagation. If moreover uq is radially symmetric and monotone 
non-increasing in |x|, then we get a clearer result: u(x,t) > 0 for all x £ R and t > 0. 


Remark 

(i) By constructed weak solution we mean that it is the limit of the approximations process that produces the 
result of Theorem 11.31 

(ii) We point out that part (a) of the theorem would still be true when m > 3 once we supply an existence 
theory based on approximations with solutions of regularized problems. 


Organization of the proofs 


• In Section [3] we derive useful energy estimates valid for all m > 1. Due to the differences in the computations, 
we will separate the cases m / 2,3 and m = 3. 

• In Section |4j [ 5 ] and |fi] we prove the existence of a weak solution of Problem 0 as the limit of a sequence of 
solutions to suitable approximate problems. The range of exponents is 1 < m < 3. 


• Section [7] deals with the property of finite speed of propagation for m > 2. See Theorem 1 7. 1| 

• In Section [8] we prove the infinite speed of propagation for m £ (1,2) in the one-dimensional case. This 
section introduces completely different tools. Indeed, we develop a theory of viscosity solutions for the integrated 
equation v t +\v x \ m ~ l (— A) 1_s w = 0, where v x = u the solution of (1.1), and we prove infinite speed of propagation 
in the usual sense for the solution v of the integrated problem. 


Though we do not get the same type of infinite propagation result for 1 < m < 2 in several spatial dimensions, 


the evidence (partial results and explicit solutions) points in that direction, see the comments in Section 10 


2 Functional setting 


We will work with the following functional spaces (see [19]). Let s £ (0,1). Let F denote the Fourier transform. 
We consider 

H s (R N ) = lu:L 2 {R N ): [ {1 + \f\ 2s )\Fu(0\ 2 df <+00 


with the norm 

l|w||ff»(R«) = IMIl2(r«) + [ l£| 2s |-7 r 'u(£)| 2 d£- 

Jr n 
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For functions u £ H S (R N ), the Fractional Laplacian is defined by 

(-A^Or) = C N> , P.V. f U{ - X) ~^ B dy = CT~\ |£| 2s (.Fu)), 

Jr n \x — y\ 

where C N>a = tt-^ 2s+n ^T(N/2 + s)/T(-s). Then 

ll w ll/f s (K jv ) = II u IIl 2 (r w ) + c||( — ^)'^ 2 ' u IIl 2 (r n )- 

For functions it that are defined on a subset Cl C with u = 0 on the boundary dCl, the fractional Laplacian 
and the H S (R N ) norm are computed by extending the function u to all R^ with u = 0 in R w \ Cl. For technical 
reasons we will only consider the case s<l/2inlV = l dimensional space. 

The inverse operator (—A) _s coincides with the Riesz potential of order 2s that will be denoted here by JC S . It 
can be represented by convolution with the Riesz kernel K s : 

JC s [u\ = K s * u, K s (x) = j |a;|^ (Ar ~ 2s) , 

c{N, s) 

where c(N, s) = n N ' 2 ~ 2s T(s)/T((N — 2s)/2). The Riesz potential KL S is a self-adjoint operator. The square root 
of K. s is JC a / 2 , i.e. the Riesz potential of order s (up to a constant). We will denote it by Hs := (/C^) 1 / 2 . Then 
Us can be represented by convolution with the kernel K s / 2 . We will write K, and % when s is fixed and known. 
We refer to [2B] for the arguments of potential theory used throughout the paper. 

The inverse fractional Laplacian K. s [u\ is well defined as an integral operator for all s £ (0,1) in dimension 
N > 2, and s £ (0,1/2] in the one-dimensional case N = 1. We extend our result to the remaining case 
s £ (1/2,1) by giving a suitable meaning to the combined operator (V/C s ). The details concerning this case will 
be given in Section |6.5[ 

For functions depending on x and t, convolution is applied for fixed t with respect to the spatial variables and 
we then write u(t) = 


2.1 Functional inequalities related to the fractional Laplacian 


We recall some functional inequalities related to the fractional Laplacian operator that we used throughout the 
paper. We refer to m for the proofs. 


Lemma 2.1 (Stroock-Varopoulos Inequality). Let 0 < s < 1, q > 1. Then 


( 2 . 1 ) 


[ M 9 2 v(—A) s vdx > ^ 

/ R« 9 


(-A) S / 2 M 9/2 


dx 


for all v £ L q (R N ) such that (— A) s v £ L q (R N ). 

Lemma 2.2 (Generalized Stroock-Varopoulos Inequality). Let 0 < s < 1. Then 


( 2 . 2 ) 


whenever ip' = ('L') 2 . 


f ip(v)(—A) s vdx > ( (—A) s/,2 , F(u) 

Jr n J r« 


dx 


Theorem 2.3 (Sobolev Inequality). Let 0<s<l (s<| if N = 1). Then 


(2.3) 




(-A ) s/2 f 


where the best constant is given in m page 31. 
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2.2 Approximation of the Inverse Fractional Laplacian (—A) s 

We consider an approximation IC e s as follows. Let K s (z ) = c/v jS |z|~( Ar ~ 2s ' ) the kernel of the Riesz potential 
1C a = (— A) _s , 0<s<l(0<s< 1/2 if N = 1). Let p e (x) = e~ N p(x/e), e > 0 a standard mollifying 
sequence, where p is positive, radially symmetric and decreasing, p £ C^°(M. N ) and f RN pdx = 1. We define the 
regularization of K s as K\ = p e * K s ■ Then 

(2.4) Kl[u\ = Kl*u 

is an approximation of the Riesz potential /C s = (—A) _s . Moreover, 1C S and K. e s are self-adjoint operators with 
K. s = ( T~L S ) 2 , K,\ = (H e s ) 2 . Also, p = a * a where a has the same properties as p. Then, we can write H\ as the 
operator with kernel K s /2 * cr e . That is: 



u K. e s [u\dx = 



\H e s [u]\ 2 dx. 


Also H\ commutes with the gradient: 


VHl[u]=Hl[Vu}. 


3 Basic estimates 


In what follows, we perform formal computations on the solution of Problem (1.11, for which we assume smooth¬ 
ness, integrability and fast decay as \x\ —> oo. The useful computations for the theory of existence and propaga¬ 
tion will be justified later by the approximation process. We fix s £ (0,1) and m > 1. Let u be the solution of 


Problem (1.1) with initial data Uq > 0. We assume u > 0 for the beginning. This property will be proved later. 


Conservation of mass: 


(3.1) 


f u(x,t)dx= [ u t dx = f V ■ (« m 1 \7IC s [u})dx = 0. 

dt Jrn J r n J r n 


• First energy estimate: The estimates here are significantly different depending on the exponent m. There¬ 
fore, we consider the cases: 

Case m = 3: 


d_ 

dt 



\ogu(x, t)dx 


[ *dx = 

J R w u 



S7u ■ V/C s [u] 



\X7H s [u]\ 2 dx. 


Therefore, by the conservation of mass (3.11 we obtain 


(3.2) 


d 

dt 


(u — log u)dx = 



\X7H s [u]\ 2 dx. 


Case m ± 3: 

d_ 

dt 


f u 3 m ( x , t)dx = (3 — m ) f u 2 m u t dx = (3 — to) f u 2 m \ 7 ( u m 1 V/C s [w])dx 

J r w Js. N J r ™ 

= — (3 —m)(2 —m) f Vu ■ S7K, s [u]dx = -C [ \\/U s [u]\ 2 dx. 

Jr n J r ™ 


Here C — (3 — m) (2 — m) is negative for m £ (2, 3) and positive otherwise. 
If ?n > 3 or 1 < m < 2 then 


d 

dt 


[ u 3 ~ m dx = -\C\ [ \\m s [u]\ 2 dx. 

J R w J R N 
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If 2 < m < 3 then 


d 

dt 


or equivalently 

• Second energy estimate: 

(3.3) 


[ u 3 ~ m dx = \C\ f \\7H s [u]\ 2 dx, 

Jr™ Jr™ 

-y- f u-u 3 ~ m dx= -\C\ f \S7'H s [u\\ 2 dx. 
dt Jr™ Jrw 


/R™ 


• L°° estimate: We prove that the L° 
Xq of u at time t = to, we have 


I \H s [u](x,t)\ 2 dx = [ 'H a [u]{'H a [u]) t dx = [ K, s [u\u t dx 
A dt J R ™ J R ™ Jrn 

/C s [u]V • (u m ~ 1 VIC s [u])dx = — [ u m ~MM[u]\ 2 dx. 

Jr™ 


norm does not increase in time. Indeed, at a point of maximum 


u t = (to — l)u m i Vu-Vp+u m i A/C s [u]. 

The first term is zero since Vu(x 0 ,fo) = 0. For the second one we have —A K, s = {— A)(—A) _s = (—A) 1_s so 
that 

A/C s [u](x 0 ,to) = -(-A) 1_s u(x 0 , t 0 ) = -c U ^^^ N ^ 2 (i-s) d d - °> 
where c = c(s, N) > 0. We conclude by the positivity of u that 

u t (xo, to) = u rn ~ 1 (x 0 ,t 0 )A/C s [u](x 0 ,t 0 ) < 0. 

• Conservation of positivity: we prove that if uq > 0 then u(t) > 0 for all times. The argument is similar to 
the one above. 


• L p estimates for 1 < p < oo. The following computations are valid for all m > 1, since p + m — 2 > 0: 

— [ u p (x,t)dx=p f u p_1 V • {u m ~ 1 '\7K. s {u])dx 
dt Jr™ Jrn 

= -p [ u m - 1 V(u p - 1 )-VIC s [u]dx = - p ( p ~ 1 \ f V(u p+m_2 ) • \7K. s [u]dx 
Jrn m 4- p — 2 J^n 

- PlP-V.f u p+m ~ 2 AIC s [u}dx = - P ^ P ~ f u p+m - 2 (-A) 1 - s u dx 
Jr™ m + p — 2 J R ™ 


< 


ui -\~ p — 2 

Mp -1) 


(m + p — l) 2 


1-3 m+p-1 

(—A) 2 u 2 


dx, 


where we applied the Stroock-Varopoulos inequality (2.11 with r = m+p+ 1. We obtain that / u p (t)dx is 

Jr™ 

non-increasing in time. Moreover, by Sobolev Inequality (2.3) applied to the function / = u ( m +P- 1 )/ 2 5 we obtain 
that 


— / u p (x, t)dx < — 
dt J R ™ 


Mp -!) 


(m+p — l) 2 5 2 _ s 


'R™ 


W(m+p-l) 

\u(x,t)\ ™- 2 + 2 » dx 


with the restriction of s > 1/2 if N = 1. 















4 Existence of smooth approximate solutions for m £ (l,oc) 


Our aim is to solve the initial-value problem (1.1) posed in Q = R™ x (0, oo) or at least Qt = R 7 '' x (0, T), with 


parameter 0 < s < 1. We will consider initial data u o £ L 1 (R JV ). We assume for technical reasons that uq is 
bounded and we also impose decay conditions as |ar| —> oo. 


4.1 Approximate problem 

We make an approach to problem 0 based on regularization, elimination of the degeneracy and reduction of 
the spatial domain. Once we have solved the approximate problems, we derive estimates that allow us to pass 
to the limit in all the steps one by one, to finally obtain the existence of a weak solution of the original problem 
0 . Specifically, for small e, 6 ,/x £ (0,1) and R > 0 we consider the following initial boundary value problem 
posed in Qt,r = {x £ B R { 0), t £ (0,T)} 


(PcS/j-R) 


(l/i)t = SAUi + V • for (x,t) £ Qt,r , 


t/i(a;,0) = u 0 {x) 
Ui (x, t) = 0 


for x £ B r { 0), 

for x £ dB R ( 0), t > 0. 


The regularization tools that we use are as follows, uq = uo,e,R is a nonnegative, smooth and bounded approxi¬ 
mation of the initial data u o such that ||wo||oo < ||uo||oo for all e > 0. For every jx > 0, c£ : [0, oo) — > [0, oo) is a 
continuous function defined by 


(4.1) 


dn(v) = {v + n) 


m —1 


The approximation of K. e s of IC S = (—A) s is made as before in Section[2j The existence of a solution U\(x, t) to 
Problem (P^rr) can be done by classical methods and the solution is smooth. See for instance [27] for similar 


arguments. 

In the weak formulation we have 

( 4 . 2 ) r 


Ui(<f>t — 6A<j))dxdt — f ( d^{U\)WK e s \Ui]S7 cjidxdt + f uo(x)(f>(x 1 0)dx = 0 

J 0 J B r J Br 

valid for smooth test functions <fi that vanish on the spatial boundary dB R and for large t. We use the notation 
Br = Br{ 0). 


Notations. The existence of a weak solution of problem (1.1 1 is done by passing to the limit step-by-step in 


the approximating problems as follows. We denote by U\ the solution of the approximating problem (P e Sfj,R I 
with parameters e, <5, /z, R. Then we will obtain 1 / 2 ( 2 ;, t) = lim e _>o Ui(x, t). Thus U 2 will solve an approximating 


problem (Ps^r) with parameters <5,/i,H. Next, we take U 3 = hin^^oo 1 / 2 ( 2 ;) that will be a solution of Problem 
(Pfj,s I, t /4 := lim^o 1 / 3 ( 2 ;, t) solving Problem (| // [) . Finally we obtain u(x,t) = lim^^o 1 / 4 ( 2 ;, t) which solves 


problem ( 1 . 1 ). 


4.2 A-priori estimates for the approximate problem 


We derive suitable a-priori estimates for the solution U\(x,t) to Problem (P t s^ R ) depending on the parameters 
e, 6 ,/x,R. 

dU 

• Decay of total mass: Since U\ > 0 and U\ = 0 in DBr, then —— < 0 and so, an easy computation gives us 

on 


d 

dt ,j B 


(4.3) 


f U±(x,t)dx = 5 f AU\ dx + I V • (d M (t/i)V/C c [l/i])d: 
Jbr Jb r Jb r 

= / f 

JdB R on JdB R on 


da < 0. 
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We conclude that 



U\{x, t)dx < 



for all t > 0. 


• Conservation of L°° bound: we prove that 0 < U\(x,t) < ||uo||oo- The argument is as in the previous 
section, using also that at a minimum point AUi > 0 and at a maximum point AU\ < 0. Also at this kind of 
points we have that 


Vd ll (U 1 ) = d' ti (U 1 )VU 1 =0. 


• Conservation of non-negativity: U\{x,t) > 0 for all t > 0, x S Br. The proof is similar to the one in the 
previous section. 


4.2.1 First energy estimate 


We choose a function F M such that 


F m ( 0) = F;(0) = 0 and F'^u) = l/d». 

Then, with these conditions one can see that F )L (z) > 0 for all z > 0. Also F^Ui) and F'^Ui) vanish on 
dB r x [0,T], therefore, after integrating by parts, we get 


(4.4) 


^ I F^Ujdx =-6 [ 

m J b r Je 


I Wr 


b r 


d^{u) 


■dx — 


\S7Ul[Ui]\ 2 dx, 


1 /2 

where 'H t =K <i . This formula implies that for all 0 < t < T we have 

(4.5) [ F^{U\{t))dx + S [ [ [7m \ dxdt + / / \^K[Ui]\ 2 dxdt = f F^u^dx. 

Jb r Jo JBr J o JBr JBr 

This implies estimates for |V"H'(17i)| 2 and 5\S7Ui\ 2 /d^(Ui) in L l {Q TR ). We show how the upper bounds for 
such norms depend on the parameters e, S , R , /i and T. 

The explicit formula for F ;J is as follows: 

1 


m +- d 3 ~ m } 


1 


(2 — m) (3 — to) 

F ^) = \ - log (1 + (Cfr/M)) + t/i/M, 

(f/r+^log (1 + (C/r/zx)) - C/r, 


2 — TO 


/x ” m 17i for to ^ 2,3, 
for to = 3, 
for to = 2. 


From formula (4.4) we obtain that the quantity F /J ,(Ui(x,t))dx is non-increasing in time: 

0< / F^(Ui(x,t))dx < / F^(u 0 )dx, Vt > 0. 

-'Br ■'Br 

Then, if we control the term f B F /x (uo)dx : we will obtain uniform estimates independent of time t > 0 for the 
quantity 


s [ [ dxdt + [ I \vw a [Ui]\ 2 dxdt. 

Jo Jb r d^,{Ui) J 0 J Br 


These estimate are different depending on the range of parameters in. 
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• Uniform bound in the case m £ (1, 2). We obtain uniform bounds in all parameters e, R, 5, p. for the energy 
estimate (4.51, that allow us to pass to the limit and obtain a solution of the original problem (1.1). By the 
Mean Value Theorem 


/ Ffj,(uo)dx < 
Jb r 


1 


< 


(2 — to) (3 — to) J b 

1 


Jb r 


2 — TO 
1 

< - 

2 — TO 

Our main estimate in the case to £ (1, 2) is: 


/ (u 0 + p) 2 m uodx 

Jb r 


(ll M o||oo + 1) 


2 — 771 


uodx. 


I«. N 


(4.6) 


where C\ = C\ ( m,uo) = 


s[ [ ^^dxdt+ [ [ \W'H\\U\\\ 2 dxdt < Ci, 
Jo Jb r dftyu i) J o J Br 


(2 — to ) 


|tto||oo + l) 2 m ||' u o||L 1 (R™)- This is a bound independent of the parameters 


e, 6, R and p. 

• Upper bound in the case m £ (2,3). 


/ F^(uo)dx = - 
Jb r 


1 


(to — 2) (3 — to) J b 


[ [(u 0 + ^- m -^- m ]dx+^—^- m f ■ 
Jb r m - a j Br 


uodx 


< 


1 


TO — 2 


uodx < 


1 


TO — 2 


uodx. 


This upper bound will allow us to obtain compactness arguments in e and R for fixed fj,. We will be able to 
control J Br F li (uo)dx — J Ba F^iUiit^dx uniformly in /x, after passing to the limit as e —► 0 and R — > oo, due to 
a exponential decay result on the solution at time t £ [0, T] that we will prove in Section [5] and the conservation 
of mass. 

Remark. These techniques do not apply in the case m > 3 because even an exponential decay on the solution 
is not enough to control the terms in the first energy estimate. 


4.2.2 Second energy estimate 


Similar computations to (3.3) yields to the following energy inequality 


— / mu^dx < -6 [ i vuimfdx-[ (u.+^r-^iciiu^dx. 
A ai Jb r Jb r Jb r 


I Br 

This implies that, for all 0 < t < T we have 

rT 


(4-7) | / \W a [Ui(t)]\ 2 dx + S j / \V-Hl[U 1 ]\ 2 dx+[ f (U, + ^) m ~ 1 \VKl[U 1 ]\ a dx < \ f |^[S 0 ]| 2 ^. 

" J Br J 0 J Br J 0 J Br " J Br 


Note that the last integral is well defined as long as uq £ T 1 (K JV ) D L° 


5 Exponential tail control in the case m > 2 

In this section and the next one we will give the proof of Theorem |1.3[ Weak solutions of the original problem 
are constructed by passing to the limit after a tail control step. We develop a comparison method with a suitable 
family of barrier functions, that in m received the name of true supersolutions. 
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Theorem 5.1. Let 0 < s < 1/2, m > 2 and let U\ be the solution of Problem (PcS^r). We assume that U\ is 
bounded 0 < Ui(x,t) < L and that uq lies below a function of the form 

V 0 (x) = Ae~ aM , A,a> 0. 


If A is large, then there is a constant C > 0 that depends only on ( N , s, a, L , A) such that for any T > 0 we will 
have the comparison 


lh(x, t ) < Ae Ct ~ a W for all x € R N , 0 < t < T. 


Proof. • Reduction. By scaling we may put a = L = 1. This is done by considering instead of U\, the function 
Ui defined as 

(5.1) Ui{x,t) = LU 1 (ax,bt), b = L m ~ 1 a 2 ~ 2s , 

which satisfies the equation 

(Ui)t = S.AU, + V.(d^(f7i)V/Cr(f/i)), 

with <5i = a 2s 6/L m ~ 1 . Note that then U\(x, 0) < A\e~^ with A\ = A/L. The corresponding bound for Ui(x,t) 
will be U\{x,t) < A/L e c ^~W with C x =C/b = C (L m ~ 1 a 2 ~ 2s )~ 1 . 

• Contact analysis. Therefore we assume that 0 < Ui(x,0) < 1 and also that 

Ui(x, 0) < Ae~ r , r =\x\> 0, 

where A > 0 is a constant that will be chosen below, say larger than 2. Given constants C, e and g > 0, we 
consider a radially symmetric candidate for the upper barrier function of the form 

U{x,t) = Ae Ct ~ r + hAe vt , 


and we take h small. Then C will be determined in terms of A to satisfy a true supersolution condition which 
is obtained by contradiction at the first point (x c , t c ) of possible contact of u and U. 

The equation satisfied by u can be written in the form 


(5.2) {U 1 ) t = 6 AU 1 + {rn-l)(u + p) m - 2 VU 1 -Vp+(U 1 +p,) m - 1 Ap, p = JC e s [U 1 }. 


We will obtain necessary conditions in order for equation (5.2) to hold at the contact point ( x c ,t c ). Then, we 
prove there is a suitable choice of parameters C, A, 77 , h, /i such that the contact can not hold. 


Estimates on u and p at the first contact point. For 0 < s < 1/2, at the first contact point (x c ,t c ) we 
have the estimates 


d r Ui = -Ae Ctc r °, AUi < Ae Ctc Ta , (Ui) t > ACe ct ° r ° + hrjAe ??tc . 


Since we assumed our solution u is bounded by 0 < u < 1, then 

(5.3) Ui(x c ,t c ) = Ae ct *- r * + hAe < 1 . 

Moreover, from m we have the following upper bounds for the pressure term at the contact point for 0 < s < 1 / 2 : 


(5.4) 


Ap(x c ,t c ) < Ki, (-d r p){x c ,t c ) < K 2 . 


Note that we are considering a regularized version of the p used in m- Of course the estimates still true (maybe 
with slightly bigger constants) since U\ is regular. 


Necessary conditions at the first contact point. Equation (5.2) at the contact point ( x c ,t c ) with r c = \x c 
implies that 


A Ce ct e - rc + hr]Ae r,tc < SAe Ct °~ r <= + (m - 1) (t/i(a; c ,t c ) + p ) m ~ 2 (-Ae Ct ^){d r p)+ 
+ (Ui{x c ,t c ) + p) m ~ 1 Ap. 
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We denote £ := r c + (r] — C)t c . Using also (5.4) with K = max/ifi, K 2 }, we obtain, after we simplify the previous 
inequality by Ae Cta ~ rc , 


C + hpe^ < <5 + (to — 1) (U\(x c , t c ) + p) m 2 K + (U\(x c ,t c ) + p) m 2 (1 + he^ + ^-e r ° Ctc )K, 
and equivalently 

C + aje^ < 5 + K (u(x c , t c ) + /z) m_2 (m + he^ + ^e rc ~ Ct °^ . 

We take C = n and ^- < h. Then 
A ~ 

C + hCe r ° < 6 + K (i U x {x c ,t c ) + p) m ~ 2 (m + he r ‘ + he r °~ Ctc ) . 


Moreover, 


By (5.3) we have that 


Since m > 2, then 


C + hCe rc <5 + /i(Ui(x c ,t c )+Ai) m_2 (m + 2/ie rc )- 
M < Ui (x c , t c ) + p < 1 + p. 


C + hCe r ° <S + K(l + p ) m ~ 2 (to + 2he r °). 

This is impossible for C large enough such that 

(5.5) C > <5 + toA'(1 + p) m ~ 2 and C > 2K (1 + p) m ~ 2 . 

Since p < 1 and S < 1, then we can choose C as constant, only depending on to and I\. 


□ 


Theorem 5.2. Let 1/2 < s < 1, m > 2. Under the assumptions of the previous theorem, the stated tail estimate 
works locally in time. The global statement must be replaced by the following: there exists an increasing function 
C(t) such that 

(5.6) u(x,t) < Ae c ^*~ a \ x \ for all x £ and all 0 < t < T. 


Proof. The proof of this result is similar to the one in m but with a technical adaptation to our model. When 
N > 2, 1/2 < s < 1, the upper bound A p{x c ,t c ) < K 0 at the first contact point holds. Moreover, in [TO], the 
following upper bound for (—d r p)(x c ,t c ) is obtained, 

(-d rP )(x c ,t c ) < k\+K^\ um\\' q \\Uim ( &- 1), \ 

where 1 < q < N/(2s — 1). We know that 11 fTi (t) | |oo < 1 and before the first contact point we have that 
U\{x,t) < Ae ct e~W, hence ||{7i(f)||i < K^Ae Ct . Therefore, if we consider K = max{Ko, KiK.\] q } we have 
that 

( 5 . 7 ) Ap(x c ,t c ) < K, (-d r p)(x c ,t c ) < K + KA 1/q e Cto/q . 

Using this estimates in the equation we obtain 

C + hge^ < 6 + K(m — 1 ) (u(x c ,t c ) + p) m ~ 2 (1 + A 1 ^ q e Ctc ^ q )+ 

+ K{u{x c , t c ) + p) m ~ 2 ( 1 + he? + j^e ra ~ Cta ). 
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We put C = T], h = p/A and use that p < u(x c , t c ) + p < 1 + p to get 


C + hCe r ° < S + KA}/ q {m - 1) (1 + p ) m ~ 2 e ct ^ q + K( 1 + p ) m ~ 2 (m + 2/ie rc ). 

We consider p < 1. The contradiction argument works as before with the big difference that we must restrict 
the time so that e Gtc ^ q < 2, which happens if 

t c <T 1 = (qlog2)/C. 


Then 

C + hCe Ta < S + 2 m ~ 1 KA 1 / q (m - 1) + 2 m ~ 2 Km + 2 m ~ 2 Khe r ‘. 

Since A > 1 and S < 1, and hence 2 m ~ 1 KA 1 ^ q (m— 1) + 2 m ~ 2 Km < 2 m ~ 1 KA 1 J q {2m — 1), we get a contradiction 
by choosing C such that: 

C = 2 m KA x l q m > 6 + 2 rn ~ 1 KA 1 ^{2m - 1). 

We have proved that there will be no contact with the barrier 

B 1 (x,t) = Ae Ct ~W 

for t <Ti = ci A~ x l q where Ci = ■ 

We can repeat the argument for another time interval by considering the problem with initial value at time 71, 
that is, 

U\{x, Xi) < Ae CTl ~W = A ie ~ lxl where A i = Ae CT \ 

and we get U(x,t) < e Clt ~for T\<t <T 2 = c\A~^ x J q e~ GTx J q where C\ = Ce GTl / q . In this way we could find 
an upper bound to a certain time for the solution depending on the initial data through A. 

When N = 1, 1/2 < s < 1, the operator d r p and A p are considered in the sense given in Section 6.5 III 


6 Existence of weak solutions for m £ (1,3) 

6.1 Limit as e —>■ 0 


We begin with the limit as e —> 0 in order to obtain a solution of the equation 
(Pshr) ( U 2 ) t = SAU 2 + V • (d^U 2 )VIC s [U 2 }). 


Let U\ be the solution of (P £ <5 m _r )■ We fix S, p and R and we argue for e close to 0. Then, by the energy formula 
(4.61 and the estimates from Section 4.2 we obtain that 

(6.1) sf [ < C(p,m,u 0 ), [ [ dxdt < C(p,m,u 0 ), 

Jo Jb r \ u i + W Jo Jb r 

valid for all e > 0. Since ||E/i||oo < ||uo||oo for all e > 0, then 

f f \VUi\ 2 dxdt < (7(/i,m,u 0 )(||uo||oo + l) m- \ Ve > 0. 

Jo Jbh 


We recall that in the case m £ (1, 2) the constant C is independent of p, that is C = C(m, uq). 

I. Convergence as e -> 0. We perform an analysis of the family of approximate solutions {U\) e in order to 
derive a compactness property in suitable functional spaces. 

• Uniform boundedness: U\ £ L°°(Qt,r)i an d th e bound |]Ui(t)]|.t,oc(R«) < | |'ito||i,°“(R N ) is independent of e, 6 ,p 
and R for all t > 0. Moreover ||/7i(t)|| L i( R w) < ||x*o|| z. 1 (k jv ) f° r all t > 0. 
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• Gradient estimates. From the energy formula (6.1) we derive 


U! £ L 2 ([0,T] : £ L 2 ([0,T] : L 2 (R N )) 

uniformly bounded for e > 0. Since VH|[?7i] is “a derivative of order 1 — s of f/i”, we conclude that 
(6.2) Ui € L 2 ([0,T],i? 1_s (M Ar )). 


Estimates on the time derivative ( U \) t : we use the equation (P^s^r ) to obtain that 


(6.3) 


(t/r)t ei 2 ([0,T] iJT" 1 ^)) 


as follows: 

(a) Since t/i G L 2 ([0,T] : Hq(Br)) we obtain that At/i € L 2 ([0,T] : H~ 1 (R N )). 

(b) As a consequence of the Second Energy Estimate and the fact that U\ £ L°°(Qt ), we have that d M (£/i)V/Cs[C/i] £ 


L 2 ([0,T] : L 2 {R N )), therefore V • JC S [U i])) £ L 2 ([0,T] : TJ" 1 ^)). 

, allow us to 

apply the compactness criteria of Simon, see Lemma |9.3| in the Appendix, in the context of 


Now, since 11(t^i)i11([ 0 ,r] : j e/- i+-(m~)) < rl/2 ||(^i)«lU;qo,Tl : g-i+«(R *))> expressions (6.2) and (6.3 


H 1 ~ s (R n ) C L 2 ^) C H~ 1 (R n ), 


and we conclude that the family of approximate solutions (U\) is relatively compact in L 2 ([0,T] : L 2 (R N )). 
Therefore, there exists a limit (U\ )e,s,n,R -> (^2)5, /x, it as e -5> 0 in L 2 ([0,T] : ^(R^)), up to subsequences. Note 
that, since (U\) e is a family of positive functions defined on Br and extended to 0 in R N \ Br, then the limit 
U 2 = 0 a.e. on R^ \ Br. We obtain that 


(6.4) 


U\ ^ U 2 in L 2 ([0,T\ : L 2 (B R )) = L 2 (B R x [0,T]). 


II. The limit is a solution of the new problem (Ps^r). More exactly, we pass to the limit as e —>• 0 in the 


definition (4.2) of a weak solution of Problem (P e s^R) and prove that the limit U 2 (x,t) of the solutions Ui(x,t ) 
is a solution of Problem (Ps^r). The convergence of the first integral in (|4.2[) is justified by (6.4) since 


’Br 


(Ui — C/ 2 ) (0* — 5A<f>)dxdt 


< \\Ul - U 2 1 |l 2 (Br x [0,T]) | |0t — SA(f>\\ L 2 (B R x[0,T])- 


Convergence of the second integral in (4.2) is consequence of the second energy estimate (4.7) as we show now. 
First we note that 

\\(Ui+fi) 2 V/CJ[f7i]|| L 2 (Bi?x ( 0jT) ) < C 

for some constant C > 0 independent of e. Then, Banach-Alaoglu ensures that there exists a subsequence such 
that 

(U 1 + aO^V/CM ^ v in L 2 {B r x (0, T)) weakly . 

Moreover, it is trivial to show that {U\ + /i) -222 ^ (U 2 + /x)” 2 ^ in L 2 (Br x (0,T)). Then 




(Ui +/x)’" 2 


(Ui + n) 2 


-VICs[Ui 


£->0 


(U 2 + n) 2 


— in L\B r x (0, T)). 


In particular we get that there exists a limit of V/C^[C/i] as e — > 0 in any L p (Br x (0, T)) with 1 < p < 00 . Now 
we need to identify this limit. The following Lemma shows that \7IC e s [Ui} V/C s [[/ 2 ] in distributions, and so 
we can conclude convergence in L 2 {Br x (0,T)). 
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Lemma 6.1. Let s £ (0,1) (0 < s < 1/2 if N = 1). Then 
(1) K*[lh\ ^ JC S [U 2 ] in L\B r x (0, T)). 


( 2 ) 


n 

to J Be 


IClp^Vi/tdxdt 


e—XO 


/7 

to JBe 


IC s [U 2 ]'S/il>dxdt for every ip £ C^{Qt)- 


Proof. For the first part of the Lemma, we split the integral as follows, 

[ [ (>C e s [U 1 ]-IC s [U 2 ])dxdt= f f (!C e a [U 1 ]-lC a [U 1 ])dxdt + f [ (K, s [Ui]-lC s [U 2 })dxdt. 

Jo J Br Jo J b r Jo J Br 

Note that 1C s [Ui\ = K S *U\ with K s £ L\ oc { R w ) and tC L s \Ui] = K e s * U\ with K f s = p e *K s where p e is a standard 
modifier. Then the first integral on the right hand side goes to zero as e —> 0. The second integral goes to zero 
with e as consequence of (6.4). 

The second part of the Lemma is just a corollary of the first part. 

r T r 


'0 J Br 


{K\[Ui\-K a [Ui\)V^dxdt 


< ||VV>||oo||/C'[L r i] -/C s [J 7 2 ]|U hBHxto.r))- 


□ 


The remaining case N 


1, s £ (1/2,1) will be explained in Section 6.5 


We conclude that, 




d^U^VtCHU^cfdxdt 



d l _ l {U 2 )ytC s \U 2 ]y(fdxdt, as e —> 0. 


Note that we can obtain also that VH e s [Ui] V"H S [J72] in L 2 (Br x (0,T)) using the same argument. This 
allows us to pass to the limit in the energy estimates. 


The conclusion of this step is that we have obtained a weak solution of the initial value problem ( Ps^r ) posed 
in Br x [0, T] with homogeneous Dirichlet boundary conditions. The regularity of U 2 , [t/ 2 ] and /C s [C/ 2 ] is as 

stated before. We also have the energy formulas 


(6.5) 


tB R 


F l _ l (U 2 (t))dx ■ 


I vt / 2 


Br d,ft 


(U2) 


dxdt - 


\S 7 H s [U 2 ]\ 2 dxdt = / F^(u 0 )dx. 


0 J Br 


' Br 


j \U s [U 2 (t)]\ 2 dx+5 f f \Wn s [U 2 ]\ 2 dxdt+ f f (f/ 2 + /z) m - 1 |V/C s [t/ 2 ]| 2 dxdt 

Br Jo J Br Jo J Br 


4 


Br 


|^ s [u 0 ]| 2 ^- 


We do not pass now to the limit as S —> 0, because we lose Ft 1 estimates for U 2 and we deal with the problem 
caused by the boundary data. Therefore, we keep the term <5A£/ 2 . 


6.2 Limit as R —>• 00 

We will now pass to the limit as R —» 00 . The estimates used in the limit on e in Section [T2] are also independent 

on R. Then the same technique may be applied here in order to pass to the limit as R —> 00 . Indeed, we get 

that U 3 = lim [/ 2 in L 2 (M. N x (0, T)) is a weak solution of the problem in the whole space 
R —>00 

(P^) {U 3 ) t = 6AU 3 + V-((U 3 + iJt) m - 1 V1C B [U 3 ]) x£R n , t> 0. 

This problem satisfies the property of conservation of mass, that we prove next. 
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Lemma 6.2. Letuo £ L 1 (M. N )r\L°° (R N ). Then the constructed non-negative solution of Problem (P^s) satisfies 

(6.6) / U 3 (x,t)dx = / uo(x)dx for all t > 0. 

J I* JR N 


Proof. Let ip £ a cutoff test function supported in the ball B 2 r and such that = 1 for \x\ < R, we 

recall the construction in the Appendix |9.2| We get 


[ (U 3 ) t tpdx = 6 f U 3 Aipdx - j ( U 3 + ■ S7<pdx = h + h- 

J B2R J B2R J B2R 

Since U 3 (t ) £ L 1 (K JV ) for any t > 0, we estimate the first integral as I\ = 0(R~ 2 ) and then I\ —> 0 as R —> 00 . 
For the second integral we have 


h = 


f IC s [U 3 ]\t ■ ((C/ 3 + ip) dx, 

J B2R 

h = (m - 1) f IC S [U 3 ](U 3 + p) m - 2 Vu ■ Vipdx + j IC S [U 3 }(U 3 + fi^A ipdx = I 21 + . 

J Bn R j Bo. Ft 


22 - 


Since VC / 3 £ L 2 (R N ) and U 3 £ L°°(R N ), 
I/21 1 <C\\{U 3 + p) m - 2 \\, 


\ 1/2 


' B 21 


\WU 3 \ 2 dxJ 


' B 2 R 


\ICs[U 3 }\ 2 \V(p\ 2 dx 


1/2 


1 


2N 


Now Vy> = V p> £ L p with p > N, so we need IC S [U 3 \ £ L q for q < 2 1 — 

— Nj2 N — ^ 

since IC[U 3 } £ L q for q > q 0 = N/(N — 2s), and q 0 < 2N/(N — 2) if 4s < N + 2. So, since p> N, 

1/9 / r \ Vp 


which is true 


I/21I < c(f \VJC s [U 3 ]\ q dx) ([ \Vp\ p dx 

\J B2R / \J B2R 


c 


For I 22 , we will use the same trick of the previous section, 


i/p 

R~ p dx ) < CR "V s 0. 


I 22 — 


[ IC S [U 3 ] [(C /3 + /i ) m_1 - /i” 1 - 1 ] Aipdx + /i m_1 f K s [U 3 ]A(pd: 

J B2R J B2R 


x = I, 


221 ~r ^ 222 - 


Now, 


I 222 — 


M m_1 f U 3 K s [Ap\dx = p m ~ 1 \\U 3 \\iO(R~ 2+2s ) 0, 

J B2R 


where we use the fact that K.A has homogeneity 2 — 2s > 0 as a differential operator. Also, 


-C 221 — 


' B2R 


f[£)U 3 K a {U 3 )A<pdx, 


where /(s) = s m 1 and £ £ [/x, /x + U 3 {x)\. Again, since U 3 £ L°°, there exists a global bound for /'(£), that is, 
/'(£) < («i- — l)max{/x m_2 , (/x + HC/alloo)” 1-2 } and so integral / 2 21 —>• 0 as i? —> 00 (details could be found in 

nni). 


In the limit R —> 00 , (p = 1 and we get (6.6). 


□ 
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Consequence. The estimates done in Section |4~2] can be improved passing to the limit R — > oo, since the 
conservation of mass (6.6) eliminates some of the integrals that presented difficulties when trying to obtain 
upper bounds independent of /i. Therefore, we compute the following terms in the energy estimate (6.5). 

For m / 2,3 we have 


[ F /1 (uo)dx — j F l _ l (U 2 )dx = 

JBr JBr 

= C [ [(w 0 + y) 3 ~ m - y 3 ~ m ]dx - 

JBr 

- C f [(U 2 + y) 3 ~ m - y 3 ~ m ]dx + 

J Br 


1 


.2 —m 


2 — TO 
1 

2 — rn 


ugdx 


'Br 


U 2 dx 


(6.7) 


C [ [(w 0 + y) 3 ~ m - y 3 ~ m }dx - C [ [(J7 3 + y) 3 ~ m ~ y 3 ~ m ]dx , 

J R« JR N 


as R —> oo. We use the notation C = 


(2—m)(3— m) ’ 


For to = 3 we have 


[ F l j - (uo)dx — [ F^(U 2 )dx = 
J Br J Br 


’Br 


Br 

Wo 


y 


1 


log 1 H- ) dx H— / woda; + / log 1 H-da;-/ L^da; 




( 6 . 8 ) 


/R« 


£7 3 


A* 


log 1 H-— ) dec — / log 1 H-I da; as i? —>■ oo. 


'Bb 

w 0 


t/ 2 




1 


A 1 as £ 




The following theorem summarizes the results obtained until now. 

Theorem 6.3. Let m > 1 and Uq € L 1 (M N ) n L°°(R N ) be non-negative. Then there exists a weak solution U 3 
of Problem (P^s I posed in x (0, T) with initial data uq. Moreover, U 3 £ L°°([0, 00 ) : L 1 (K Ar )), and for all 
t > 0 we have 


U 3 (x,t)dx = / uo(x)dx 

J R« 


and ||J7 3 (.,t)|| 00 < 11Wo11 00 - The following energy estimates also hold: 

(i) First energy estimate: 

• If to = 3, 

-dxdt + f f |V7A s [b r 3 ]| 2 da:dA 
’ Jo J R™ 

< 


(6.9) 


4 ' |VC/ 3 | 2 


0 V {U3 + y) 2 


log (— + 1^ dx 

\T J 


f log (?«*> +iU. 

Jr n \ y J 


• If 2,3 and 

( 6 . 10 ) 


|W 3 | 2 


/ (TT J- \m-l dxdt + f f \^'Hs[U 3 ]\ 2 dxdt + 

0 JR. N (JJ 3 + y) J 0 J r n 

+c f [{U 3 (t) + y) 3 ~ m - y 3 ~ m ] dx < cf [(w 0 + y) 3 ~ m - y 3 ~ m ] 

J r» Jm N 


dx 


where C = C{m) = (2 _ m) 1 (3 _ m) ■ 
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(ii) Second energy estimate: 


1 

2 



\n s [u 3 (T)}\ 2 dx+s [ 

Jo 



\S7H s [U 3 ]\ 2 dxdt + / 

JO J R 


R N 


1 

< - 
~ 2 


iHsluaj^dx. 


6.3 Limit as /r —>• 0 


Similarly to the previous limits we can prove that t /4 
will be a solution of problem 

(■ Ps) ( U A ) t = 6 AU 4 + V • (t/ 4 m - 


= lim U 3 in 

fi—tO 




L 2 ^ x 


x G 




(0, T)) when m G (1,3). Then t / 4 

t > 0. 


In order to pass to the limit, we need to find uniform bounds on g > 0 for terms 3 and 4 of the energy estimates 
(6.9) and ( 6 . 10 ). 

Uniform upper bounds 

• Case m G (1, 2). By the Mean Value Theorem, 


(to — 2)(3 — 


Jtx -T / [(«o + M) 3 ” m - t?- m ] dx < 

)(3 - to) J r n 

z (||uo||oo + l) 2 - m 


1 


TO — 2 


/ 


(to — 2) 

uodx. 


[ (u 0 + g ) 2 m u 0 dx 

Jr n 


This bound is independent of g. 

• Case to G (2,3). The function /(£) = <^ 3_m is concave and so f(U 3 + g) < f(g) + f(U 3 ). In this way, 


1 


(2 — to) (3 — to) 


/ [(U 3 (t) + g) 3 ~ m - g 3 ~ m ] dx < 

Jr n 


1 


(2 — m)( 3 — m) 


[ U 3 (t) 3 ~ m dx. 

J R« 


The last integral is finite due to the exponential decay for U 3 that we proved in Section [5j In this way, the last 
estimate is uniform in g. 

The limit is a solution of the new problem ( Pg I. The argument from Section 6.1 does not apply for the 
limit 

(6.11) [ [ (U 3 +y) m - 1 VIC s [U 3 }V(t>dxdt'^ [ [ £/J* -1 V/C a [If 4 ]V0dsdt. 

Jo Js. N Jo J R w 


In order to show that this convergence holds, we note that from the first energy estimate we get that 


VH,[U 3 ] Gi 2 ((0,T) :T 2 (R jv )) 

uniformly on g. Then V/C. s [U 3 \ = 'H S [\7JJ S [U 3 ]\ G L 2 ((0,T) : H S (R N )). Since for any bounded domain f!, H S (U) 
is compactly embedded in L 2 (fl) then X7IC S [U 3 ] —>• \7K, S \U^\ as g —> 0 in L 2 (fl). Then we have the convergence 
since U 3 G L°°(R N ) and (j) is compactly supported. 

Remarks. • In the case m = 2 the corresponding term is f RJV U 3 log - (t /3 + g)dx which is uniformly bounded if 
U 3 has an exponential tail. This has been proved by Caffarelli and Vazquez in m- We do not repeat the proof 
here. 

• The case m > 3 is more difficult since we can not find uniform estimates in g > 0 for the energy estimates 
that allow us to pass to the limit. 


( 6.11 
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6.4 Limit as 5 —>■ 0 


We will prove that there exists a limit u = lim^o U 4 in L 2 (R N x (0, T)) and that u{x,t) is a weak solution to 


Problem (1.1). Thus, we conclude the proof of Theorem 1.2 stated in the introduction of this chapter. 


We comment on the differences that appear in this case. From the first energy estimate we have that 


/R« U. 


-dxdt < C(m,uo), 


which gives us that 5 \ 7 U 4 G L 2 (Qt) since U 4 G L°°(Qt)- Then, as in Section 6.1 we have that SAU 4 G 7F _1 (K W ) 
uniformly in 5. Also V(L r ™~ 1 V/C, i [t/ 4 ]) G H~ 1 (R N ) as before. Then ( U 4 ) t G iJ” 1 (K Ar ) independently on S. 
Therefore we use the compactness argument of Simon to obtain that there exists a limit 

Ui(x,t) -+ u(x,t) L 2 ((0,T) x R n ). 


Now we show that u is the weak solution of Problem (1.1). It is trivial that S f n J RN U 4 Acf> — >■ 0 as 6 — > 0. On 
the other hand, V/C s [[/ 4 ] = ’H S [V'H S [U 4 ]} G L 2 oc (Qt) uniformly on S > 0 since V'H./Ki] G L 2 (Qt) uniformly on 
S > 0. In this way, V/C s [I/ 4 ] has a weak limit in L 2 oc {Qt)- As in Lemma 6.1 (2) we can identify this limit and 
so on, V/C s [t/ 4 ] —► V/C s [it] weakly in L 2 oc (Qt) as <5 —>■ 0 and therefore 


U™~- L VJC s [U 4 ]V(j)dxdt 


( S ->0 


since C/™ 1 


Jo Js. N 
r m_1 in Lf oc (Q T ) as <5 -S> 0. 


/0 Jr n 


u m - L VJC s [u]X7(j)dxdt. 


6.5 Dealing with the case N — 1 and 1/2 < s < 1 

As we have commented before, the operator KL S is not well defined when N = 1 and 1/2 < s < 1 since the kernel 


11—2 s 


licp does not decay at infinity, indeed it grows. It makes no sense to think of equation (6.12) in terms of a 


pressure as before. This is maybe not very convenient, but it is not an essential problem, since equation 
can be considered in the following sense: 

(6.12) u t {t 7 x) = V • (u ro - 1 (V/C B )[u]) for x G R N , t > 0, 

where the combined operator (V/C s ) is defined as the convolution operator 

(V/C s )[tt] := (VAT S ) * u with K s (x ) = 


11 —2s ' 


Other authors that dealt with N = 1 have considered operator (V/C s ) before. They use the notation V 2s 1 to 
refer to it. Note that 


VK a ( x) = (-1 + 2 s)c s 


x 


,|3-2s ’ 


and so, VA' S G L} oc ( 


for N = 1 and 1/2 < s < 1. Moreover, (V/C s ) is an integral operator in this range. As 


in Subsection 2.2 the operator (V/C s ) is approximated by (V/C s ) e defined as 

(V/C s ) e [u] = (VI< s y * u where (VA s ) e = p e * (VA S ). 

Note that 

(6.13) (VA s ) E ^VIv s in Lj oc (R), 


since VA S G L/ oc (R). It is still true that 


(V/C s )[u] = U s [VAsM], 
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since the operator TL S is well defined for any s £ (0,1) even in dimension N = 1. 

In this way, almos t all the arguments from Section [6] apply by replacing V(/C s (u)) for (V/C s )(u). The only 
exception is Lemma 6.1 where the weak L 2 (M) limit of X7K s \Ui] is identified. This argument is replaced by the 
following Lemma: 


Lemma 6.4. Let N=1 and 1/2 < s < 1. Then 


' 0 J Br 


Ui(VJCsY[ ip] dxdt 




' 0 J Br 


U 2 (VJC s M dxdt \/if G C™{Qt). 


Proof. 


/o JB h 


t/i(V/C s ) e [i/] - U 2 {VIC s )[ip] dxdt = 


10 J B, 
r T r 


10 JBp 


(Ui — f7 2 )(V/C s ) e [V’] dxdt 

u 2 ((vic s y[ii>}-(vjc s M)dxdt. 


The first integral on the right hand side goes to zero with e since 11 (\7/C s ) e ['0] | |z.°°(k) < K for some positive 
constant K which does not depend on e and U\ —> U 2 as e —» 0 in L 2 (Br x (0,T)). The second integral also 
goes to zero as consequence of (6.13) and the fact that U 2 £ L°°(R) uniformly on e. □ 


7 Finite propagation property for m £ [2,3) 


In this section we will prove that compactly supported initial data uq{x) determine the solutions u(x,t) that 
have the same property for all positive times. 


Theorem 7.1. Let m > 2. Assume u is a bounded solution, 0 < u < L, of Equation (1.1) with 1C = (—A) _s with 
0<s<l (0 < s < 1/2 if N = 1), as constructed in Theorem 1.3 Assume that uq has compact support. Then 
u(-,t) is compactly supported for all t > 0. More precisely, if 0 < s < 1/2 and uq is below the ”parabola-like” 
function 

U 0 {x) = a{\x\ - b) 2 , 

for some a,b > 0, with support in the ball Bb( 0), then there is a constant C large enough, such that 


u(x, t) < a(Ct — (|.t| — b )) 2 . 


Actually, we can take C(L, a) = C (1,1) L m 2+s a 2 s . For 1/2 < s < 1 a similar conclusion is true, but C = C(t) 
is an increasing function oft and we do not obtain a scaling dependence of L and a. 


Proof. The method is similar to the tail control section. We assume u(x, t) > 0 has bounded initial data 
uq(x) = u(x, t 0 ) < L, and also that uq is below the parabola Uq(x) = a(|a;| — b ) 2 , a, b > 0. Moreover the support 
of Uq is the ball of radius b and the graphs of Ug and Uq are strictly separated in that ball. We take as comparison 
function U(x, t) = a(Ct — (|a;| — b )) 2 and argue at the first point in space and time where u{x, t) touches U from 
below. The fact that such a first contact point happens for t > 0 and x y oo is justified by regularization, as 
before. We put r = |a;|. 


By scaling we may put a = L = 1. We denote by ( x Cl t c ) this contact point where we have u(x c , t c ) = U(x c , t c ) = 
(b + Ct c — |^c|) 2 - The contact can not be at the vanishing point |a:/(t c )| := b + Ct c of the barrier and this will 
be proved in Lemma 7.2 We consider that x c lies at a distance h > 0 from |a;/(t c )| = b + Ct c (the boundary of 
the support of the parabola U{x,t) at time t c ), that is 


b + Ct c — |*c| = h > 0. 
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Note that since u < 1 we must have \h\ < 1. Assuming that u is also C 2 smooth, since we deal with a first 
contact point (x c , t c ), we have that u = U, X7(u — U) = 0, A (u — U) < 0, (it — U) t > 0, that is 

u(x c ,t c ) = h 2 , u r = — 2h, Ait < 2AT, u t >2Ch. 

For p = /C s (it) and using the equation it t = (m — l)it m-2 Vw • Vp + it m_1 Ap, we get the inequality 


(7.1) 


2Ch < 2 (to — l)/i 2m 3 ( — p r + — Ap ) , 


where p r and Ap are the values of p r and Ap at the point (x c ,t c ). In order to get a contradiction, we will use 
estimates for the values of p7 and Ap already proved in [10] (see Theorem 5.1. of [10]) 

(7.2) -P7< K ± + K 2 h 1+2s + K 3 h, Ap < I< 4 . 

Therefore, inequality (7.1| combined with the estimates fl7 implies that 

(7.3) 2 C < 2(m - 1 )h 2m ~ 4 (R\ + K 2 h 1+2s + Kh ) , 

which is impossible for C large (independent of h), since m > 2 and \h\ < 1. Therefore, there cannot be a 
contact point with h ^ 0. In this way we get a minimal constant C = C(N, s) for which such contact does not 
take place. 


Remark: For m < 2, we do not obtain a contradiction in the estimate (7.3), since the term K\h 2rn 4 can be 
very large for small values of \h\. 

• Reduction. Dependence on L and a. The equation is invariant under the scaling 
(7.4) u(x,t) = Au(Bx,Tt) 

with parameters A,B,T> 0 such that T = A m ~ l B 2 ~ 2s . 

Step I. We prove that if u has height 0 < u(x,t) < 1 and initially satisfies u(x,0) = uq(x) < (|x| — b ) 2 then 
u(x, t) < U(x, t) = (Ct — (|x| — b )) 2 for all t > 0. 

Step II. We search for parameters A, B , T for which the function u is defined by (7.41 satisfies 

0 < u(x, t) < L , u(x, 0) < a(|a:| — b) 2 . 

An easy computation gives us 

A = L, AB 2 = a, b = b/B. 

Moreover, by the relation between A, B and T we obtain A = L, B = (a/L) 1 / 2 and then T = L m ~ 2 +s a 1 ~ s . 
Then u(x, t) is below the upper barrier U(x, t) = a(Ct — (|x| — b )) 2 where the new speed is given by 

C = CA^B 1 - 28 = CL m ~i +s a!“ s . 


• Case 1/2 < s < 1. The proof relies on estimating the term d r p at a possible contact point. This is independent 
on m and it was done in [TDl. G 


Lemma 7.2. Under the assumptions of Theorem 7.1 there is no contact between u(x, t) and the parabola U ( x , t), 
in the sense that strict separation of u and U holds for all t > 0 if C is large enough. 

Proof. We want to eliminate the possible contact of the supports at the lower part of the parabola, that is the 
minimum \x\ = Ct + b. Instead of analyzing the possible contact point, we proceed by a change in the test 
function that we replace by 


U e (x,t) = 


(Ct — (|cc| — b )) 2 + e(l + Dt) for \x\ < b + Ct , 


-(1 + Dt), 


for \x\ > b + Ct. 
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The function U £ is constructed from the parabola U by a vertical translation e(l + Dt) and a lower truncation 
with 1 + Dt outside the ball {|x| < b + Ct}. Here 0 < e < 1 is a small constant and D > 0 will be suitable 
chosen. 

We assume that the solution u(x,t) starts as it(x, 0) = uq(x) and touches for the first time the parabola U £ at 
t =t c and spatial coordinate x c . The contact point can not be a ball {|x| < b + Ct} since U £ is a parabola here 
and this case was eliminated in the previous Theorem |7.1| Consider now the case when the first contact point 
between u(x,t ) and U £ (x,t) is when |x c | > b + Ct c . At the contact point we have that u = U £ , V(it — U £ ) = 0, 
A (it — U £ ) < 0, (it — U e ) t > 0. In this region the spatial derivatives of XJ £ are zero, hence the equation gives us 

De = {e{l + Dt c )) m ~ 1 &p, 

where A p is the value of A p = (—A) 1-S w at the point (x c ,t c ). Since e is small we get that the bound u(x,t) < 
U±(x,t) is true for all |x| < R N . This allows us to prove that that A p is bounded by a constant K. We obtain 
that De < (e(l + Dt c )) m ~ 1 K. Since m > 2 and e < 1, this implies that 

D < (1 + Dt c ) m ~ 1 K. 

We obtain a contradiction for large D 1 for example D = 2K , and for 

tc<Tc= ^K ( 2l/(m_1) - X ) ■ 

Therefore, we proved that a contact point between u and U £ is not posible for t < T c , and thus u(x,t ) < U £ (x,t) 
for t < T c . The estimate on t c is uniform in e and we obtain in the limit e —> 0 that 

u(x,t ) < U(x,t) = (Ct — (|x| — b )) for t < /^ 2 1 A m ^ 1 ) _ x'j _ 

As a consequence, the support of u(x,t) is bounded by the line |x| = Ct + b in the time interval [0,T C ). The 
comparison for all times can be proved with an iteration process in time. 

• Regularity requirements. Using the smooth solutions of the approximate equations, the previous conclusions 
hold for any constructed weak solution. 

□ 

Remark. The following result about the free boundary is valid only for s < 1/2 and for solutions with bounded 
and compactly supported initial data. The result is a direct consequence of the parabolic barrier study done 
in the previous section. Since that barrier does not depend explicitly on m if m > 2, the proof presented in 
m is valid here. By free boundary J-B(u ) we mean, the topological boundary of the support of the solution 
5(it) = {(x,t) : u{x,t) > 0}. 

Corolary 7.3 (Growth estimates of the support). Let uq be bounded with uq(x) = 0 for |x| > R for some 
R > 0. If (x, t) £ J-B{u) then x < R + Ct 1 /( 2 -‘ 2s ) ; where C = C'(||tio|| 0 o! s )- 


7.1 Persistence of positivity 

This property is also interesting in the sense that avoids the possibility of degeneracy points for the solutions. 
In particular, assuming that the solutions are continuous, it implies the non-shrinking of the support. Due to 
the nonlocal character of the operator, the following theorem can be proved only for a certain class of solutions. 

Lemma 7.4. Let u be a weak solution as constructed in Theorem \ 1. 3| and assume that the initial data uq(x) is 
radially symmetric and non-increasing in |x|. Then u(x,t) is also radially symmetric and non-increasing in |x|. 


Proof. The operators in the approximate problem ( P e Sfj,R^ are invariant under rotation in the space variable. 
Since the solution of problem (P £ Siir) is unique, then we obtain that u(x,t) is radially symmetric. □ 
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Theorem 7.5. Let u be a weak solution as constructed in Theorem m and assume that it is a radial function 
of the space variable u(\x\,t) and is non-increasing in |a"|. If Uq(x) is positive in a neighborhood of a point Xq, 
then u{x o, t ) is positive for all times t > 0. 

Proof. A similar technique as the one presented in the tail analysis is used for this proof, but with what we 
call true subsolutions. Assume uq(x) > c > 0 in a ball Br(x o). By translation and scaling we can also assume 
c = R = 1 and :ro = 0. Again, we will study a possible first contact point with a barrier that shrinks quickly in 
time, like 

(7.5) U(x,t)=e- at F(\x\), 

with F : R>o —> R>o to be chosen later and a > 0 large enough. Choose F(0) = 1/2, F(r) = 0 for r > 1/2 
and F'(r) < 0 for all r £ R> 0 . The contact point (x c ,t c ) is sought in B 1 / 2 ( 0 ) x (0, 00 ). By approximation we 
can assume that u is positive everywhere so there are no contact points at the parabolic border. At the possible 
contact point ( x c ,t c ) we have 

^'(‘^'C; t c ) — tl(x c ,t c ), Ut(x c: tc) A U t (x c ,t c ) — Ctl/(x c ,t c ), 

Vu(i c ,t c ) = Vl/(i c , t c ) = e~ atc F'(\x c \)e r , e r = x c /\x c \. 

We recall the equation 

u t = (m- 1 )u m - 2 VuVp + it™' 1 Ap. 

Then at the contact point ( x c ,t c ) we have 

-aU = U t > ut = (m- 1 )U m ~ 2 VUVp + U^Ap, 

where A p = A p(x c ,t c ). Then 

—ae _a * c T(|a; c |) > (m - l)e-^ m - 2 ^F{\x\) m - 2 e- a ^F'(\x c \)pf + 

According to m we know that the term F'(\x\) p r > 0 and Ap is bounded uniformly. Therefore 

—ae~ atc F(\x c \) > 

Simplifying and using that m > 2, A p is bounded uniformly and also F is bounded, we obtain 

a < — e _a ( m_2 ) tc F(|x c |) m_2 Ap < AT" a(m - 2)tc < K. 

This is not true if a > K and we arrive at a contradiction. O 


Remark. There exist counterexamples on the persistence of positivity property when the hypothesis of Theorem 


data with not connected support. 


7.5 are not satisfied. In [101 (Theorem 6.2) the authors construct an explicit counterexample by taking an initial 


8 Infinite propagation speed in the case 1 < m < 2 and N = 1 

In this section we will consider model <0 

(8.1) d t u = d x ■ (■ u m ~ 1 d x p ), p = (—A )~ s u, 

for x £ R, t > 0 and s £ (0,1). We take compactly supported initial data u 0 > 0 such that u 0 £ Li oc m- We 
want to prove infinite speed of propagation of the positivity set for this problem. This is not easy, hence we 
introduce the integrated solution v, given by 

(8.2) v(x,t) = f u(y,t)dy> 0 for t > 0, x £ R. 

J — OO 
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Therefore v x = u and v(x, t ) is a solution of the equation 

(8.3) d t v = ~\v x \ m - 1 (-A) a v , 


in some sense that we will make precise. The exponents a and s are related by a = 1 — s. The technique of 
the integrated solution has been extensively used in the standard Laplacian case to relate the porous medium 
equation with its integrated version, which is the p-Laplacian equation, always in ID, with interesting results, 
see e. g. [25] . The use of this tool in [2j for fractional Laplacians in the case m = 2 was novel and very fruitful. 
We consider equation (8.3) with initial data 


(8.4) 


v(x, 0) = vq(x) := / uq(x) dx for all x £ 


Note that v(x, t) is a non-decreasing function in the space variable x. Moreover, since u(x, t) enjoys the property 
of conservation of mass, then v{x,t) satisfies (see Figure [T]) 

lim v(x,t) = 0, lim v(x,t) = M 

x—>—oo x —>-+oo 


for all t > 0. We devote a separate study to the solution v of the integrated problem (8.3) in Section 8.3 The 


validity of the maximum principle for equation (8.3) allows to prove a clean propagation theorem for v. 


Theorem 8.1 (Infinite speed of propagation). Let v be the solution of Problem (8.3)-(8.4), and assume 
that uq > 0 is compactly supported. Then 0 < v(x, t) < M for all t > 0 and x £ R. 


The use of the integrated function is what forces us to work in one space dimension. The result continues 
the theory of the porous medium equation with potential pressure, by proving that model (8.1) has different 
propagation properties depending on the exponent m by the ranges m > 2 and 1 < m < 2. Such a behaviour 
is well known to be typical for the classical Porous Medium Equation ut = A u m , recovered formally for s = 0, 
which has finite propagation for m > 1 and infinite propagation for m < 1. Therefore, our result is unexpected, 
since it shows that for the fractional diffusion model the separation between finite and infinite propagation is 
moved to m = 2. 


Proof of Theorem 1.4 PART b). This weaker result follows immediately. In fact, in Theorem 8.1 
that v(x,t) defined by (8.21 is positive for every t > 0 if x £ 
arbitrary far from the origin such that u(x, t ) > 0. 


we prove 

Therefore for every t > 0 there exist points x 


If moreover, uq is radially symmetric and non-increasing in |a;| and u inherits the symmetry and monotonicity 
properties of the initial data as proved in Lemma |7.4[ This ensures that u can not take zero values for any x £ K. 
and t > 0. 


□ 


8.1 Study of the integrated problem 


• Connection between Model (|8.1[) and Model 


We explain how the properties of the Model (8.1) with N = 1 can be obtained via a study of the properties 


of the integrated equation (8.3). We consider equation (8.1) with compactly supported initial data uq such 


that uq > 0. Let us say that supp uq C [—/?.,!?], where R > 0. Therefore, the corresponding initial data to be 
considered for the integrated problem is vq(x) = f* Uo(y)dy, for all i£l. Then vq : R —> [0, oo) and has the 
properties 


(8.5) vq{x) = 0 for x < —R , vq(x) = M for x > R, v' 0 (x) > 0 for x £ (—R,R), 

where K > 0 is fixed from the beginning and M = f R uo(x)dx is the total mass. 
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Figure 1: Typical compactly supported initial data for models (8.1) and (8.3). 


8.2 Regularity 


Proposition 8.2. The solution v :[0,T]xl-> [0, oo) of Problem (8.3) defined by formula 

[X 

v(x,t) = / u(y,t)dy is continuous in space and time. 

J — OO 


Proof. I. Preliminary estimates. Since v x (x, t) 
the estimates of Section |6.1| we have the following: 


u{x,t) 1 where u is the solution of Problem (1.1), then by 


• v x = u £ L°°([0,T] : L°°(R)), therefore v £ L°°([0,T] : Lip(R)), where Lip(R) is the space of Lipschitz 
continuous functions on R. In particular, v £ L°°([0, T] : Lip(£? B )) for every Br C R. 


• We have (v t ) x = u t = d x (u m ~ 1 d x (—A)~ s u) in the sense of distributions. Then v t £ L 2 ([0,T] : L 2 (B)) for 
every set BcR, with \B\ < +oo. The proof is as follows. The first equality holds in the distributions sense, 
that is 


f f Vt.<p x dxdt:=— f ( v(<p x )tdxdt= f f u m 1 d x {— A) s uip x dxdt, Vip £ C^°(R x [0, T]). 
Jo Jr Jo Jr Jo Jr 

This implies that vt = u m ~ 1 d x (~A)~ s u a.e. in R. Then, using the second energy estimate (3.3), we obtain 

ll v *lli 2 ([0,T]:L 2 (S)) = ll^” 1 A) 

r T r 


< u 


i m— 1 


v \d x {—A) s u\~dxdt < +oo. 


'o Jb 


II. Continuity in time. Let (xo,to) £ R x [0, T]. Let (x,ti) £ R x [0, T] and h := x — xq. Let B = [xo,Xi]. 
Then 

\v{x 0 ,h) - v{x 0 ,t 0 )\ < |u(a;o,ti) - v(x,h)\ + \v(x,t 0 ) - v(aj 0 ,i 0 )| + \v(x,h) - v(x,t 0 )\. 
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We know v £ Lip a .(K); let L the corresponding Lipschitz constant. Then 


1 f x 

■j(xo,h)-v(x 0 ,t 0 )\<2Lh+ T / \v{y,h) - v{y,t 0 )\dy 

h Jx „ 


^ 2Lh — 


Vtdt 


dy < 2 Lh 2 


\v t \ dydt 


' Xq Jto 


< 2 Lh + -|-B| 1 *”|f 1 — ^o| 1/,2 ||wt|lz,2([o,T]:L 2 (s) 


= 2 Lh 


|*i - *o| 1/2 

ft 1 / 2 


ll' U t|lz, 2 ([0,T]:L 2 (B)- 


Optimizing, we choose ft ~ ^ tl . —, that is ft ~ (t± — *o) 3 ^ 2 , and we obtain that 


\v(x 0 ,h) - v(x 0 ,t 0 )\ < K\ti - * 0 | 1/3 - 

This estimate holds uniformly in x £ R and it proves that v(x,t) is Holder continuous in time. In particular 
v £ C'([0,T] : <7(R)). □ 


8.3 Viscosity solutions 

Notion of solution. We define the notions of viscosity sub-solution, super-solution and solution in the sense 
of Crandall-Lions m The definition will be adapted to our problem by considering the time dependency and 
also the nonlocal character of the Fractional Laplacian operator. For a presentation of the theory of viscosity 
solutions to more general integro-differential equations we refer to Barles and Irnbert [lj. 

It will be useful to make the notations: 

USC(Q) = {upper semi-continuous functions u : Q —► R}, 

LSC(Q) = {lower semi-continuous functions u : Q —> M}, 

C(Q) = {continuous functions u : Q —> R}. 

Definition 8.3. Let v £ USC(M. x (0, oo)) (resp. v £ LSC(R x (0,oo)) ). We say that v is a viscosity sub¬ 
solution (resp. super-solution) of equation (8.3) onlx (0, oo) if for any point (xo,to) with to > 0 and any 
t £ (0, to) and any test function tp £ C 2 )!! x (0, oo)) D L°°(R x (0, oo)) such that v — tp attains a global maximum 
(minimum) at the point (cco,to) on 

Q t = R x (t 0 - r, t 0 ] 


we have that 


d t p(x 0 ,t 0 ) + \<Px{x 0 ,to)\ m 1 ((-A)“^(-,t 0 ))(a;o) < 0 (> 0). 


Since equation (8.3) is invariant under translation, the test function p in the above definition can be taken such 


that p touches v from above in the sub-solution case, resp. p touches v from below in the super-solution case. 


We say that v is a viscosity sub-solution (resp. super-solution) of the initial-value problem (8.3)-(8.4) on 
R x (0, oo) if it satisfies moreover at t = 0 

v(x,0) < limsup v(y,t ) (resp. v(x,0)> liminf v(y,t)). 

y^fx, t—► o t - > ° 

We say that v £ C( R x (0, oo)) is a viscosity solution if v is a viscosity sub-solution and a viscosity super¬ 
solution on R x (0, oo). 


Proposition 8.4 (Existence of viscosity solutions). Let u be a weak solution for Problem (1.1). Then v defined 

r _ _ 

by formula v{x 1 t) = / u(y,t)dy is a viscosity solution for Problem (8.3)-(8.4). 
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Proof. By Proposition 8.2 we know that v £ C([0,T] : C(R)). The idea is to obtain a viscosity solution by the 

f x 

approximation process. Let v e defined by v e (x,t ) = / u e (y,t)dy , where u e is the approximation of u as in 

J — OO 

Section [4] Then v e is a classical solution, in particular a viscosity solution, to the problem 

(i> e )t = SA(v e ) + |(v e ) a: | T "- 1 (-A) 1 — 

Since u e —X u, then we get that v e —> v as e — > 0 (and similarly with respect to the other parameters). The final 
argument is to prove that a limit of viscosity solutions is a viscosity solution of Problem (8.3)-(8.4). 

□ 


The standard comparison principle for viscosity solutions holds true. We refer to Imbert, Monneau and Rouy 
[23] where they treat the case m = 2 and a = 1/2. Also, we mention Jakobsen and Karlsen [24] for the elliptic 
case. 


Proposition 8.5 (Comparison Principle). Let m £ (1, 2), a £ (0, 1), N = 1. Let w be a sub-solution and W be 
a super-solution in the viscosity sense of equation (8.31. If w(x, 0) < Vq < W(x, 0), then w < W in 1 x (0, oo). 


We give now our extended version of parabolic comparison principle, which represents an important instrument 
when using barrier methods. This type of result is motivated by the nonlocal character of the problem and the 
construction of lower barriers in a desired region O Cl possibly unbounded. This determines the parabolic 
boundary of a domain of the form LI x [0, T] to be (R \ LI) x [0,T]URx {0}, where Ll C R. A similar parabolic 
comparison has been proved in [Bj and has been used for instance in mm- 


Proposition 8.6. Let m > 1 , a £ (0,1). Let v be a viscosity solution of Problem (8.3)-(8.4). Let<& : 
R such that <f> € C 2 (Ll x (0, T)). Assume that 


lx [0, oo) 


• + |4> a; | m-1 (—A)“4> < 0 forx£Ll,t£ [0,T]; 

• 4>(a;, 0) < v(x, 0) for all x £ R (comparison at initial time); 

• &(x,t) < v(x,t ) for all x £M.\Ll and t £ (0, T) (comparison on the parabolic boundary). 


Then 4>(a :,t) < v(x,t) for all x £R, t £ (0,T). 


Proof. The proof relies on the study of the difference $-»:Rx [0, oo) —X R. At the initial time t = 0we have 
by hypothesis that $(£, 0) — v(x, 0) < 0 for all x £ R. 

Now, we argue by contradiction. We assume that the function <f> — v has a first contact point (x c ,t c ) where 
x c £ Ll and t c £ (0, T). That is, ($ — v)(x c ,t c ) = 0 and (<£> — v)(x,t) < 0 for all 0 < t < t c , x £ R, by regularity 
assumptions. Therefore, (<f> — v) has a global maximum point at ( x c ,t c ) on R x (0,f c ]. Therefore, v — attains 
a global minimum at ( x c ,t c ). 

Since v is a viscosity solution and $ is an admissible test function then by definition 

®t(x c ,tc) + \^x{x c ,tc)\ rn ~ 1 {-A) a ^(x c ,t c ) > 0 , 

which is a contradiction since this value is negative by hypothesis. O 


8.4 Self-Similar Solutions. Formal approach 

Self-similar solutions are the key tool in describing the asymptotic behaviour of the solution to certain parabolic 
problems. We perform here a formal computation of a type of self-similar solution to equation (8.3), being 
motivated by the construction of suitable lower barriers. 
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Let m £ (1, 2) and a £ (0,1). We search for self-similar solutions to equation (8.3) of the form 


U(x,t) = $(\y\t ) 


which solve equation (8.31 intx (0, oo). After a formal computation, it follows that the exponent b > 0 is given 
by b = l/(m — 1 + 2a) and the profile function $ is a solution of the equation 

We deduce that any possible behaviour of the form <L(y) = c |?/|~ 7 with 7 > 1 is given by 

2a + m 


( 8 . 6 ) 


7 = 


2 — m 


The value of the self-similarity exponent will be used in the next section for the construction of a lower barrier. 
A further analysis of self-similar solutions is beyond the purpose of this paper and can be the subject of a new 
work. We mention that in the case m = 2, the profile function d> has been computed explicitly by Biler, Karch 
and Monneau in [3]. 


8.5 Construction of the lower barrier 


In this section we present a class of sub-solutions of equation (8.3) which represent an important tool in the 


proof of the infinite speed of propagation. For a suitable choice of parameters this type of sub-solution will 

give us a lower bound for v in the corresponding domain. This motivates us to refer to this function as a lower 

barrier. We mention that a similar lower barrier has been constructed in ED- 

to + 2a 1 , 

Let 7 = —- and 0 = --— — be the exponents deduced m Section 


2 — TO 


to — 1 + 2a 


8.4 


We fix Xq < 0. In the sequel we will use as an important tool a function G : 
constants C\ > 0 and C 2 > 0, we have that 

• (Gl) G is compactly supported in the interval (—Xo,oo); 

• (G2) G(x) < Ci for all i£l; 

• (G3) (—A ) s G(x) < — C 2 \x\~ ( ~ 1+2s ' > for all x < xq. 

This technical result will be proven in Lemma 9.1 of Section [9] (Appendix). 


such that, given any two 


Lemma 8.7 (Lower Barrier). Let xq < 0, e > 0 and £ > 0. Also, let G be a function with the properties 
(G1),(G2) and (G3). We consider the barrier 

(8.7) $ e (*, t) = (t + r) bl ((|cc| + £)“ 7 + G(x)) - e, t > 0, x £ M. 

Then for a suitable choice of the parameter C 2 > 0, the function 4> e satisfies 

( 8 . 8 ) ($ 6 )t + |($ e h| m_ 1 (-A) Q $ e <0 for x < x 0 , t > 0. 

Moreover, C\ is a free parameter and C 2 = C 2 (N, to, a, t). 


Proof. We start by checking under which conditions <I> e satisfies ( 8 . 8 ), that is, $ is a classical sub-solution of 


equation (8.3) in Q. To this aim, we have that 


($ e ) t +K$ e ) x r- i (-A)“$ £ = 6 7 


(t + T ) bl ~ 

(M + 0' 


(+ _L T )bj(m-1) 




(M + fl’ 


(1*1 +^)G+1)(™-1) 


(|*| +^)(7+l)(m~l) 

x| + 0- 7 ] + (-A)“g). 
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Now, by Lemma 9.2 we get the estimate (—A)“ ((|a:| + £) -7 ) < C 2 ,\x\~ ( ' 1+2a ' > for all |x| > |cco|, with positive 
constant C 3 = C 3 (N,m,a). At this step, we choose the parameter C 2 in the assumption (G2) to be at least 
C 2 > C 3 . The precise choice will be deduced later. Since 7 = (7 + l)(m — 1) + 1 + 2a, we continue as follows: 


($ e ) t+ K$ e ) x r- i (-A)“$ e 

i +7 - 


(t + r) 6 ^ 


-(C' 3 -C' 2 )|cr|-( 1+2 “) 


(|x|+0 7 ' ' (M+0 (7+1)(m_1) 

= (|*| +e)- (7+1)(m - 1} - 

• (lry(t + r) 67_1 (|a:| + e)- (1+2a) + + r) b ^ m (C 3 - C 2 )\x\~( 1+2a ^ 

< (|a;| + ^) _( 7 + 1 )(m_ 1 ) |a;| _(1+2 “ ) (try(t + r ) 67 " 1 + 7 m " 1 (t + r) b 7 m (C 3 - C 2 )) 
which is negative for all ( x,t ) £ Q , if we ensure that C 2 is such that: 

(8.9) C 2 > C 3 + 6 7 2 - m r b7(1 - m) - 1 . 

This choice of C 2 is independent on the parameters £, e. 


□ 


From now on, we will take r = 1, which will be enough for our purpose. We can now prove the main result for 
the model (8.3) which in particular implies the infinite speed of propagation of model (1.1) for 1 < m < 2 in 
dimension N = 1. 


8.6 Proof Theorem 18.II 


Let *0 < 0 fixed. We prove that v(x, t) 
with properties (8.5), satisfies 


u ior an 1 /> 


and x < xq. By scaling arguments, the initial data 


( 8 . 10 ) 


v 0 (x) > H x0 {x) = 


0 , x < Xq , 
1 , X > Xq- 


We will prove that v(x, t) > d> e (j:, t) in the parabolic domain Qt = {x < xq, t £ [0, T]} by using as an essential 
tool the Parabolic Comparison Principle established in Proposition | 8 . 6 | We describe the proof in the graphics 
below, where the barrier function is represented, for simplicity, without the modification caused by the function 
G(-) (Figure [3]). 

To this aim we check the required conditions in order to apply the above mentioned comparison result. 


Comparison on the parabolic boundary. This will be done in two steps. 


(a) Comparison at the initial time. The initial data (8.10) naturally impose the following conditions on 
At time t = 0 we have ^(zq, 0) < 0, which holds only if £ satisfies 


(8.11) £ > x 0 + e t . 

Therefore >F e (xo,0) < v 0 (x 0 ) since v 0 (x 0 ) > 0. 

(b) Comparison on the lateral boundary. Let k\ := min{r(a:,t) : x > Xo, 0 < t < Tj with k\ > 0. This 
results follows from the continuity v £ C([0,T] : C(R)) since uo^o) = 1- We impose the condition 


$ £ (a :,t) < v(x,t) for all x > xq, t £ [0, T]. 


It is sufficient to have 

(r+i) b7 (r 7 + ci) < a*. 
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The maximum value of T for which this inequality holds is 


( 8 . 12 ) 


T < 


ki 

^ + C i 


1/&7 


- 1 . 


We need to impose a compatibility condition on the parameters in order to have T > 0, that is: 
(8.13) £> (Aq-Ci) - d. 

The remaining parameter C\ from assumption (G2) is chosen here such that: C\ < k\. 


By Proposition 8.6 we obtain the desired comparison 

v(x, t) > < h e (x,t) for all ( x,t ) £ Qt- 


• Infinite speed of propagation. Let x\ < £o and t\ £ (0,T) where T is given by (8.12). We prove there 
exists a suitable choice of £ and e such that $ e (xi,ti) > 0. This is equivalent to impose the following upper 
bound on £: 


(8.14) 


£ < x\ + (fy + 1) 


b— £ 


We need to check now if there exists e > 0 such that condition (8.14| is compatible with conditions (8.11) and 
(8.13). For the compatibility of conditions (8.11) and (8.13|) we have 


£ 0 + e i < £ < xi + (ti + l) b e "i, 


that is, 
(8.15) 


e < 


(ti + 1) & - 1 


i "t 


Xo - Xi 


For conditions (8.13) and (8.14) we need 

(fci - Ci)“^ ^ < X\ + (ti + l) b e - - 


which is equivalent to 


(8.16) 


e < 


(h + i) b 


(h - Ci) -y — X\ 


Both upper bounds (8.15) and (8.16) make sense since 0 > xq > X\ and ki > C\. 


Summary. The proof was performed in a constructive manner and we summarize it as follows: Ci < k\, T 


4>(fi,:ri) > 0. 

This proofs that v(ti,xi) > 0 for any t £ (0,T). 


given by (8.12). Then by taking e the minimum of (8.15)-(8.16), £ satisfying (8.11 )-(8.13)-(8.14) we obtain that 


D 


Remark. The parameter £ of the barrier depends on e by (8.11) and (8.14| and therefore £ —>■ oo as e —► 0. 
Therefore ^(a;, t) —> 0 as e —> 0 for every (x, t ) £ Qt and we can not derive a lower parabolic estimate for v(x, t) 
in Qt- 
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Figure 2: Comparison with the barrier at time t = 0 



Figure 3: Comparison with the barrier at time t > 0 


9 Appendix 

9.1 Estimating the Fractional Laplacian 

In this section we are interested in estimating the fractional Laplacian of given functions. We recall the definition 
of the Fractional Laplacian operator 


(-A ) s u(x) = a NiS P.V. 


u{x) - u{y) 
I® — y\ N+2s ’ 


0 < s < 1, 


where a s a normalization constant given by 


N-\-2s ^ 

aN ’ s = n N/ 2r(-7V/2)' 

First, given the expression of the fractional Laplacian, we construct a function with the desired properties. 

Lemma 9.1. Given two arbitrary constants C\,C 2 > 0 there exists a function G : R —► [0, +oo) with the 
following properties: 

1. G is compactly supported. 
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2. G( x) < Ci for all x £ R. 

3. (— A) s G(x) < —G 2 |a;+ 1+2s ) for all x £ R. with d{x, supp{G)) > 1. 


Proof. Let R an arbitrary positive number to be chosen later. We consider a smooth function Gi : R —»• [0, +oo) 
such that Gi(j;) < C\ for all x £ R and supported in the interval [—1,1]. 

We define Gr(x) = G\{x/R). Therefore IIGrIIl^r) = -R||Gi||x,i(H.), Gr < C\ and G is supported in the interval 
[-R, 1?]. Then for |x| > R + 1 we have that 


(-A ) s G r {x) =<t s f 
Jr 


Gr(x ) - G R (y) 


\x-y | 


1+25 


dy = -<7« 


r 

l-R F 


Gr( V ) 


-y\ 


1+2 s 


dy 


< —a, 




< -cr,2- (1+2s) ||G 


fill+bRlF 


— (l+2s) _ _ 


= -<x s 2-( 1 + 2 s +||G 1 | 


i-(1+25) 


(j 2 l+ 2 s 

It is enough to choose R > — - to get (—A ) s Gr(x) < G 2 |x| _ ( 1+2 + Note that R implicitly depends 

fJ s||Gl||il(R) 

on Gi since ||Gi||ii(®) < 2 Gi. 

n 


Secondly, we need to estimate the fractional Laplacian of a negative power function. The following result is 
similar to one proven by Bonforte and Vazquez in Lemma 2.1 from [5] with the main difference that our function 
is G 2 away from the origin. We make a brief adaptation of their proof to our situation. 

Lemma 9.2. Let ip : R -+ (0, oo), <p = (|x| + £) -7 , where 7 > 1 and £ > 0. Then, for all |.t| > +o| > 1, we have 
that 

( 9 + K—A)V(+I < 

with positive constant G > 0 that depends only on 7 , £, s. 


Proof. Let us first estimate the L 1 norm of tp. 

/ <p{x)dx = / ip(x)dx + / ip(x)dx < / fC 1 dx + / x~' y dx 

J R J\x\<l J |#|>1 J\x\<l J\x\~>l 

r°° 9 

<2r 7 + 2 / r _7 dr = 2£ -7 +-- < G, G = G( 7 ,£). 

J 1 7 - 1 


Following the ideas of 0 Lemma 2.1, the computation of the (—A ) s ip(x) is based on estimating the integrals on 
the regions 

Ri = {y. \y\> 3N/2}, R 2 = \y-. M < \y\ < \B^(x), 

R 3 = {y: \x-y\<\x\/2}, R A = {y. \y\<\x\/2}. 


Therefore 


/ = 


(—A)V(z) = / 

Jr 

each of t 

<p(x) - <p(y) 


ip(x) - ip{y) 


t R1UR2UR3UR4 1 ^ y I 

We proceed with the estimate of each of the four integrals: 

dr 

1 1 , 1 + 2 s -dySu d p{x) - = 

’ \y\>S\x\/2 F y\ + J3\x\/2 


1+2 s 


dy. 


Ki 


j»l+ 2 s |^| 7 + 2 s 


Ki = Ad ( 7 , s). 
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11 = 


ip(x)-ip(y) <p(x) 


ir 2 


- y I 


1+2 s 


dy < 


Ko 


( 1 * 1 / 2 ) 


l+2s 


/• 3|*|/2 

l\x 1/2 dr= N" +2 ' 


K 2 = K 2 {'Y,s). 


Ill = 


f v{x) - +y) ^ ^ I, //,, /' i , 

Jr 3 \x-y\ 1+2s y ~ h I|z,00(b i*i/ 2(x)) 7 | x _ j / |<| x | /2 | a ; - y | 2s_1 27 

Ah f ^ 172 1 


|a ;++ 2 


r»2s— 1 — 


dr < 


K 3 

\x\i+ 2s 


I<3 = K 3 { 7, s). 


For IF we take use that when |y| < \x\/2 then \y — x\ > |x|/2 and |y| < |x| which implies y>(j/) > ip(x). We have 


7F < 


'\y\< 1*1/2 

Ah 


|y>(*) - +d)l. 

I® - 2 /| 


l+2s 


( 2 \ r / 2 \ 


— | T |i+ 2 «> Ah(7, s,£). 


Since 7 > 1, we can conclude that 

|(-A)V(*)| < \I\ + \II\ + \III\ + \IV\ = K 5 \x\-*- 2s + A4+- 1 - 2 * < K 6 \x r 1 ” 28 , V|x| > |cc 0 | > 1. 


□ 


9.2 Reminder on cut-off functions 


We remind the construction of cut-off functions. Let 


/(*) 


e- 1 ^, x > 0 , 
0 , x < 0 . 


Then / G C°°(M). Let 


F(x) 


/(*) 

f(x)+f{l-x)’ 


iel. 


Then F(x) =0 for x < 0, F(x) = 1 for x > 1 and F(x) G (0,1) for x G (0,1). We construct now the cut-off 
function ip : R N -7 ([ 0 , 1 ] by: 

ip(x) = F (2 — |x|), x G R^. 


Then G C'°°(K Ar ), ip(x) = 1 for \x\ < 1, y>(x) = 0 for \x\ > 2 and + 2 :) G (0,1) for \x\ G (1,2). The cut-off 
function for Br is obtained by 

<Pr{%) = 1>{x/R). 


Thus ipn G C , 00 (K JV ), <Pr{x) = 1 for |aj| < R , pr(x) = 0 for \x\ > 2 R and P>r(x) G (0, 1) for |rc| G (R,2R). Also, 
we have that V(< pr) = 0(R~ 1 ), A (<pr) = 0(R~ 2 ). 


9.3 Compact sets in the space L p (0,T]B) 

Necessary and sufficient conditions of convergence in the spaces L p (0, T; B) are given by Simon in |28j. We recall 
now their applications to evolution problems. We consider the spaces X C B C Y with compact embedding 
X B. 

Lemma 9.3. Let F be a bounded family of functions in L p (0,T-,X), where 1 < p < 00 and 3F/dt = {df/dt : 
f G J 7 } be bounded in A 1 (0, T;Y). Then the family F is relatively compact in A p (0, T; B). 

Lemma 9.4. Let F be a bounded family of functions in A°°(0,T;X) and dF/dt be bounded in L r (0,T;F), 
where r > 1. Then the family F is relatively compact in C(0,T; B). 
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10 Comments and open problems 


• Case to > 3. In this range of exponents the first energy estimate does not hold anymore. Therefore, we lose 
the compactness result needed to pass to the limit in the approximations to obtain a weak solution of the original 
problem. The second energy estimate is still true and it gives us partial results for compactness. In our opinion 
a suitable tool to replace the first energy estimate would be proving the decay of some L p norm. In that case 
we will also need a Stroock-Varoupolous type inequality for some approximation C e s of the fractional Laplacian. 
The technique of regularizing the kernel by convolution that we have used through this paper does not allow us 
to prove such kind of inequality. The idea is however to use a different approximation of the pressure term that 
is well suited to the Stroock-Varoupolous type inequalities. Let us mention jl4j where this kind of inequalities 
are proved for a wider class of nonlocal operators including Cf s . The technical details are involved and the new 
approximation may have an interest, hence we think it deserves a separate study. 


we proved a transformation 
and the fractional porous medium 


• Infinite propagation in higher dimensions for self similar solutions. In 

formula between self-similar solutions of the model (1.11 with 1 < to < 2 
equation u± + (—A ) s u m = 0. This way we obtain infinite propagation for self similar solutions of the form 
U{x,t) = t- a F{\x\t~ a > N ) in KA. This is a partial confirmation that the property of the infinite speed of 
propagation holds in higher dimensions for every solution of (1.11 with 1 < to < 2 . 


• Explicit solutions. Y. Huang reports m the explicit expression of the Barenblatt solution for the special 
value of m , m ex = (N 6s — 2 )/(N + 2s). The profile is given by 


F M (y) = X(R 2 + \y\ 2 )- {N+2s)/2 , 


where the two constants A and R are determined by the total mass M of the solution and the parameter (3. 
Note that for s = 1/2 we have m ex = 1, and the solution corresponds to the linear case, Ut = (—A) 1 ^ 2 u, 
F 1 / 2 (r) = C(a 2 + r 2 )“( Ar+1 )/ 2 . 


• Different generalizations of model (CV) are worth studying: 

(i) Changing-sign solutions for the problem d t u = V • (|it|Vp), p = (—A)" 


it. 


(ii) Starting from the Problem (CV), an alternative is to consider the problem 

t > 0, 


= V • (|m|V(—A) - s (| 


u\ m ~ 2 u)), 


X G 


t,N 


with to > 1. This problem has been studied by Biler, Imbert and Karch in [3]. They construct explicit compactly 
supported self-similar solutions which generalize the Barenblatt profiles of the PME. In a later work by Imbert 
|22|, finite speed of propagation is proved for general solutions. 

(iii) We should consider combining the above models into d t u = V(|u| m_1 Vp), p = (—A )~ s u. 

When s = 0 and m = 2 we obtain the signed porous medium equation dtu = A(|u| TO-1 it). 
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